Modeling Chandra X-ray observations of 
Galaxy Clusters using Cosmological 
Simulations 



Daisuke Nagai 1 , Andrey V. Kravtsov 2 , and Alexey Vikhlinin 3 ' 4 

1 Theoretical Astrophysics, California Institute of Technology, Mail Code 130-33, 
Pasadena, CA 91125 daisuke@caltech.edu 

2 Department of Astronomy and Astrophysics, KICP, & EFI, The University of 
Chicago, 5640 South Ellis Ave., Chicago, IL 60637 

3 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, 
MA 02138 

4 Space Research Institute, 8432 Profsojuznaya St., GSP-7, Moscow 117997, 
Russia 

1 Abstract 

X-ray observations of galaxy clusters potentially provide powerful cosmolog- 
ical probes if systematics due to our incomplete knowledge of the intraclus- 
ter medium (ICM) physics are understood and controlled. In this paper, we 
study the effects of galaxy formation on the properties of the ICM and X- 
ray observable-mass relations using high-resolution self-consistent cosmolog- 
ical simulations of galaxy clusters and comparing their results with recent 
Chandra X-ray observations. We show that despite complexities of their for- 
mation and uncertainties in their modeling, clusters of galaxies both in ob- 
servations and numerical simulations are remarkably regular outside of their 
cores, which holds great promise for their use as cosmological probes. 



2 Testing X-Ray Measurements of Galaxy Clusters 

X-ray observations with Chandra and XMM-Newton enable us to study prop- 
erties of the ICM with unprecedented detail and accuracy and provide impor- 
tant handles on the ICM modeling and associated systematics. Their superb 
spatial resolution and sensitivity enable accurate X-ray brightness and tem- 
perature measurements at a large fraction of the cluster virial radius and also 
make it simple to detect most of the small-scale X-ray clumps. Despite this re- 
cent observational progress, the biases in the determination of the key cluster 
properties remain relatively uncertain. 
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Fig. 1. Mock Chandra images of one of the relaxed (left) and unrelaxed (right) 
simulated clusters at z=0. The detectable extended X-ray sources, indicated by 
ellipses, are detected and masked out from the analysis. The physical size of the 
images is 5 ft" 1 Mpc. 

We therefore assess the accuracy of the X-ray measurements of galaxy 
cluster properties using mock Chandra analyses of cosmological cluster simu- 
lations and analyzing them using a model and procedure essentially identical 
to that used in real data analysis [1, 2]. The comparison of the true and de- 
rived cluster properties provides an assessment of biases introduced by the 
X-ray analysis. We examine the bias in mass measurements separately for 
dynamically relaxed and non-relaxed clusters, identified based on the overall 
structural morphology of their Chandra images, which mimics the procedure 
used by observers. The typical examples of systems classified as relaxed or 
unrelaxed are shown in Figure 1. To check for any redshift dependence in 
such biases, we also analyze the simulation outputs at z — and 0.6. The 
simulations and analysis procedures are fully described in [3]. 

Figure 2 illustrates that the X-ray analysis provides accurate reconstruc- 
tion of the 3D properties of the ICM for the nearby, relaxed clusters. The 
strongest biases we find are those in the hydrostatic mass estimates, which 
is biased at a level of about 13% at r 50 o c even in the relaxed clusters. We 
find that the biases are primarily due to additional pressure support provided 
by subsonic bulk motions in the ICM, ubiquitous in our simulations even in 
relaxed systems [4, 5, 6]. These biases are related to physics explicitly missing 
from the hydrostatic method (e.g., turbulence), and not to deficiencies of the 
X-ray analysis. Gas fraction determinations are therefore biased high. The 
bias increases toward cluster outskirts and depends sensitively on its dynam- 
ical state, but we do not observe significant trends of the bias with cluster 
mass or redshift. We also compute a X-ray spectral temperature (Tx), a value 
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Fig. 2. The mock Chandra analyses of one of the relaxed clusters at z = with 
A/500C = 1.41 x 10 14 /i _1 Mq. In the left and middle panels, the best-fit model (dashed 
lines) recovers well both the projected profiles and the actual 3D gas profiles. In 
the upper-right panel, the reconstructed Af gas profile (dot-dashed line) is accurate 
to a few percent in the entire radial range shown. The hydrostatic M to t estimate 
(dotted line), on the other hand, is biased low by about 5%-10% in the radial range, 
[0.2, l.Ojrsooc- The lower-right shows that measured cumulative / gas is biased high 
by ~ 10% in the radial range of [0.2, 1.0]r 5 oo c for this cluster, and it is primarily due 
to the bias in the hydrostatic mass estimate. 

derived from a single-temperature fit to the integrated cluster spectrum ex- 
cluding the core (< 0.15r 50 o c ) and detectable small-scale clumps [7, 8]. 

3 Effects of Galaxy Formation on the ICM Profiles 

Next, we investigate the effects of galaxy formation on the ICM properties and 
compare the results of simulations with recent Chandra X-ray observations of 
nearby relaxed clusters. The impact of galaxy formation on the properties of 
ICM are investigated by comparing simulations performed with and without 
the processes of galaxy formation (e.g., gas cooling, star formation, stellar 
feedback and metal enrichment), which we refer to as the cooling+SF (CSF) 
and non-radiative runs, respectively. Fig. 3 shows that gas cooling and star 
formation modify both the normalization and the shapes of the ICM profiles. 
What happens is that the removal of low-entropy gas in the inner region 
raises the level of entropy and lowers the gas density [9]. The effects are 
strongly radial dependent and increase toward the inner regions down to about 
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Fig. 3. The ICM gas density and entropy profiles of simulated clusters and Chandra 
X-ray observations of nearby (z ~ 0) relaxed clusters. 

~ O.lrsooc, inside which the observed properties are not well reproduced in 
the simulations. On the other hand, the ICM properties outside the cores 
in the cooling CSF simulations and observations agree quite well. At r§oo Cl 
both the ICM density and entropy profiles of different mass systems converge, 
indicating that the clusters are self-similar in the outskirts. Note that the non- 
radiative simulations predict overall shape of the density and entropy profiles 
inconsistent with observations. 

4 X-ray observable-mass relations 

For cosmological application, it is important to understand the relations be- 
tween X-ray observables and cluster mass. In Fig. 4, we present recent com- 
parisons of two X-ray proxies for the cluster mass — the spectral temperature, 
Tx, and the new proxy, Yx, defined as a simple product of Tx and M g [10]. 
Analogously to the integrated Sunyaev-Zel'dovich flux, Yx is related to the 
total thermal energy of the ICM. 

The M500 — Tx relation has a ~ 20% scatter in M^qq around the mean 
relation, most of which is due to unrelaxed clusters. The unrelaxed clusters 
also have temperatures biased low for a given mass. This is likely because 
during mergers, the mass of the system has already increased but only a 
fraction of the kinetic energy of merging systems is converted into the thermal 
energy of gas, due to incomplete relaxation [11]. The slope and evolution of 
the M500 — Tx relation are also quite close to the self-similar model. 

The M500 — Yx relation shows the scatter of only rj 7%. Note that this 
value of scatter includes clusters at both low and high-redshifts and both 
relaxed and unrelaxed systems. In fact, the scatter in M 500 — Yx for relaxed 
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Fig. 4. Correlation between the total mass, M$ooc and X-ray spectral temperature, 
Tx (left panel) and the integrated X-ray pressure, Yx (right panel) . Separate symbols 
indicate relaxed and unrelaxed simulated clusters, and also z=0 and 0.6 samples. The 
figures include points corresponding to three projections of each cluster. The dot- 
dashed lines are the power law relation with the self-similar slope fit for the sample 
of relaxed clusters. The dotted lines indicate 20% and 8% scatter, respectively. The 
data points with errorbars are Chandra measurements of nearby relaxed clusters. 



and unrelaxed systems is indistinguishable within the errors. Yx is therefore 
a robust mass indicator with remarkably low scatter in A/500 for fixed Yx, 
regardless of whether the clusters are relaxed or not. The redshift evolution of 
the Yx — M500 relation is also close to the simple self-similar prediction, which 
makes this indicator a very attractive observable for studies of cluster mass 
function with X-ray selected samples, because it indicates that the redshift 
evolution can be parameterized using a simple, well-motivated function. 

Finally, the results of the simulations are compared to the observational 
results. In both relations, the observed clusters show a tight correlation with 
a slope close to the self-similar value. The normalization for our simulated 
sample agrees with the Chandra measurements to ~ 10 — 15%. This is a con- 
siderable improvement given that significant disagreement existed just several 
years ago [12, 13]. The residual systematic difference in the normalization is 
likely caused by non-thermal pressure support from bulk gas motions, which 
is unaccounted for in X-ray hydrostatic mass estimates. The much improved 
agreement of simulations and observations in these relations gives us confi- 
dence that the clusters formed in modern simulations are sufficiently realistic 
and thus can be meaningfully used for interpretation of observations. The ex- 
istence of tight relations of X-ray observables, such as Yx, and total cluster 
mass and the simple redshift evolution of these relations hold promise for the 
use of clusters as cosmological probes. 
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